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The Jaynes-Cummings model (JCM) describes the interaction of a two-level atom with a single 
quantised field mode in an optical cavity. In this BSc project we derive the basic Rabi model 
Hamiltonian and show how it leads to the JCM Hamiltonian. The results are then used to analyse 
the general time evolution of a bipartite atom-field system and its properties, such as entanglement. 
We extend a scheme developed by Browne and Plenio 1^ to entangle two cavities that interact with 
an atom. Our scheme allows to produce states that are arbitrarily close to a Bell state with a success 
probability greater than 0.75, whereas the success rate of the original method tends to for perfect 
fidelity. Further, we extend both schemes by allowing the atom or the cavity to be in a thermal state 
and find that the fidelity of the produced states is reduced with increasing temperature. Finally, we 
apply our scheme to quantum teleportation, where we find that the fidelity of the teleported state 
compared to the original state also decreases as the temperature is increased. 



I. INTRODUCTION 



The Jaynes-Cummings model (JCM) is of particu- 
lar theoretical interest since it provides the simplest 
fully quantum mechanical description of the interac- 
tion between a two-level atom and a quantised field 
O [3]. In many cases, the full evolution of such a bi- 
partite system can be solved analytically. Hence the 
JCM represents a unique tool for analytically studying 
many properties of these systems that can otherwise 
only be evaluated numerically. Of particular interest 
is the entanglement that can be generated between 
the atom and the field, or even between individual 
field modes that interact with the atom, as shall be 
shown later. 

In this project we focus on deriving the JCM and ap- 
plying it to study the time evolution of the bipartite 
system, starting from various initial conditions. In 
section II we will derive the Rabi model Hamiltonian 
and show how, under the rotating wave approxima- 
tion (RWA), it leads to the JCM Hamiltonian. We 
will use effective Hamiltonian theory |3] to justify the 
RWA and show how the Hamiltonian behaves in the 
near resonant and far off-resonant limits. Section HI 
will utilise the found Hamiltonians to develop evolu- 
tion operators that describe the full dynamics of the 
system under the interaction. We also look at some 
particular examples with simple initial conditions that 
give a general feeling for the behaviour of the system 
as it evolves. 

In Section IV we make use of all the tools we devel- 
oped to study a particular system that was initially 
proposed by Browne and Plenio [T], which consists of 
an atom that is sent through two optical cavities and 
briefly interacts with each cavity mode, thereby entan- 
gling the two cavities. We propose an improvement on 
the original scheme and look at the effect of thermal 
interaction between the system and the environment, 
by allowing either the atom or the cavities to be in 
thermal states. Finally, in section V we evaluate if 
the generated state could be used for quantum tele- 
portation and how the fidelity of the teleported state 
depends on temperature. 



II. DERIVATION OF THE HAMILTONIAN 

Starting from the Hamiltonian of an electron which 
is bound to a nucleus under an external field [5], 

Hir,t) = -^[P + eA(r,t)]2-e$(r,t) + 1/(r), (1) 
2m 

where A(r, t) and $(r, t) are the vector and scalar po- 
tentials of the external field and V{r) is the Coulomb 
potential that binds the electron to the nucleus, we 
used gauge transformations, particularly the Coulomb 
gauge, to show that the Hamiltonian can be rewritten 
as 

iJ'(r,i) =i?o(r,t)-d-E(t), (2) 

where i7o(r,t) is the Hamiltonian ^ in the absence 
of an external field, d — — er is the dipole operator, 
and E(t) is the operator representing the time varying 
electric field [Ml] ■ The interaction between atom and 
field is represented by the term Hint = — d • E(i). By 
considering the matrix elements of this operator in the 
atomic basis |e) and \g), the excited and ground state 
respectively, one can show that (setting h = 1, as shall 
be done throughout the rest of this discussion) 

^»«t = A((7+ +(7_)(a-f a^), (3) 

where A is the coupling strength, depending on the 
dipole moment and the electric field, cr+ = |e) {g\ and 
= |g) (e| are the atomic transition operators, and 
the field raising and lowering operators and a have 
their usual meaning. 

Hq can be found by considering the energy of the 
atom and the field individually. By taking the zero 
of the atomic energy midway between the two lev- 
els, the atomic part of the Hamiltonian is simply 
Ha = 5 Woe's I where ojq is the energy difference be- 
tween the atomic levels and = \e) (e| — 1^) {g\ 
is the atomic inversion operator 5]. The Hamilto- 
nian for the field takes an equally simple form of 
Hp = uj{a^a +5), with w being the frequency of the 
radiation [5 . 

Combining the individual terms and dropping the 
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vacuum energy contribution in the field Hamiltonian 
(since it does not contribute to the dynamics), gives 
the Rabi model Hamiltonian 

H = -ojQ&j + Lxja) a + A((7_|_ + (7_)(a + a)). (4) 

Under certain conditions the interaction term can be 
further simplified. Using the approach described by 
James and Jerke [4] for deriving an effective interac- 
tion Hamiltonian, we showed that near resonance, i.e. 
for small detuning = ujq — io <^ ui ^ the interaction 
Hamiltonian 8J is given by 



Hi = A'^+\{(7+a + a^a) + 0{-)). (5) 

Hence if we limit our focus to systems for which 
w A,A we can neglect the O(^) term in ([s]), hav- 
ing justified the rotating wave approximation that is 
often made in the literature [21 [3] , to arrive at the 
full JCM Hamiltonian for the near resonant case, 
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III. TIME EVOLUTION 

A. Background 

The most general (pure) initial states of the atom 
and the field are given by linear superpositions of all 
the accessible states 



\m)a = Cg\g)+C,\e) 

oo 

\m)f = E^"!'^)' 



(9) 
(10) 



where the Cx are arbitrary coefficients, satisfying 
the normalisation conditions |Cgp + jCgp = 1 and 
l^nP ~ 1- "^^^ total initial state of the bipar- 
tite atom-field system is then given by 

\m)^\m)a®\m)f (n) 



= (w + A)y + wa^a + A((7+a + (7_a^). (6) 

A similar treatment for the case of large detuning, 
where w 3> A 3> A, leads to the interaction picture 
Hamiltonian 



Hr = 



A2 



(7) 



which gives the full JCM Hamiltonian for the off- 
resonant case, 

H = {uj + A)^+ oja^fa + ^ (aga^a -f- |e) (e|) . (8) 

In the next section we shall make use of these results 
to develop a framework for the time evolution of the 
atom-field system. 



1. Resonant Interaction, A = 

Allowing the C coefficients to vary with time, we di- 
rectly solved the interaction picture Schrodinger equa- 
tion 



d\m) 

dt 



Hi \m) 



(12) 



for the case of no detuning (A = 0), where the inter- 
action Hamiltonian takes the simple form 



Hi = A((7-)_a + a^aJ) 



(13) 



as can be directly seen by substituting A = into 
([5|. We found the analytical solution to (12) at an 
arbitrary time t as 



CO 

|i/,(t)) = \^(CeCn cos(Aty^rTT) - iCgCn+1 sin( Ai^ATTT)) |e) 

+ (^—iCeCn-i sin(Aiv^) + CgCn cos{Xt^/n)j \g) (g) \n) 



n=0 



(14) 



An equivalent but more versatile and powerful ap- 
proach (especially when working with density ma- 
trices, see below) is to explicitly calculate the time- 
evolution operator 



U = e 



-iH,t 



(15) 



By expanding the exponential into its power series 
and calculating the matrix elements, again using the 
Hamiltonian given in (13), we showed that the full 



evolution operator of the JCM at resonance, written 



in the atomic basis {|e) , is given by [8l ITO] 



U{t) 



cos{XtV aa'^) —ia 



sin(At V dt d) 



sin(AivddT ) 
y/ ddt 



cos{XtV a^a) 



Applying (16t to the initial state ( pT| ) 

\m) = m \m) 



(16) 



(17) 



exactly reproduces the result obtained in (14), as ex- 
pected. 
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Considering equation ( 14 ) , it is straightforward to 



analyse the types of interaction that are taking place 
in the resonant case. The first line of (14) shows that 
the excited atomic state |e) can interact with the field 
number state \n) by exciting the field into the higher 
number state \n + 1) while the atom itself drops into 
its ground state \g) to conserve the total excitation 
number. Alternatively this process can take place in 
reverse, whereby the atom gets excited \g) — > |e) by 
absorbing a photon from the field |n) — |n — 1). 



2. Off-Resonant Interaction, A ^ A 

The off-resonant interaction behaves distinctly dif- 
ferent in comparison to this. Starting again with the 
definition in \lh\ and using the off-resonant interac- 
tion Hamiltonian Q it is straightforward to find the 
corresponding evolution operator. The general initial 
state (11) is an eigenstate of this Hi. Hence the ex- 



ponential involving Hj that appears in (17) can be 



replaced by the equivalent exponential with the oper- 
ator replaced by its eigenvalues. This gives the state 
of the system in the off-resonant case at time t as 



oo 

\m) = E ^» [Cee-^'^t^+i)* |e) + C.e''^"* \g) 



where we have defined 



A 



(18) 



(19) 



for convenience. 

It is immediately obvious that ( 18 ) represents a very 



different interaction compared to the resonant case 
(14). Here, no energy or excitation exchange between 
the atom or the field take place. Only the relative 
phases are affected and rotate in time. 
In the following subsection we will use the results 
(14) and (18 1 to look at some illustrative exemplary 



systems. 



B. Two Simple Examples 

Having derived the general time evolution of an ar- 
bitrary system in the two distinct cases, we can now 
focus on some particular examples and analyse their 
properties. The coherent state \a) for which the C„ 
coefficients take the form 



Cn = e- 



(20) 



where a is related to the fields mean photon number n 
by n = |ap, is of particular interest. These states are 
in some sense the "most classical" field states, since 
they obey the classical equations of motion of the sim- 
ple harmonic oscillator and they can also serve as an 



<cr3> 

a) 




FIG. 1: a) Atomic inversion of an initially excited atom in- 
teracting with a coherent field with n = 10, plotted against 
the scaled interaction time At. b) Von Neumann entropy 
of the system. Series in this and the following figures were 
evaluated up to Cioo unless otherwise stated. 



idealisation of the light emitted by a laser [S]. 

It will also prove useful later on to define the notion of 

the density operator p{t), which for a pure state \il'{t)) 

is simply given by the outer product of the state with 

itself 



p{t)^\m)m)\^ 



(21) 



In terms of the density operator, the expectation value 
of an observable Ai of subsystem i (here either atom 
or field) is given by 



{A,) = Tr[p,Ai 



(22) 



where pi is the density operator of the subsystem, 
given by the partial trace of p over the other sub- 
system [To] , 

Equipped with these tools we can now look at an 
atom, initially in the excited state (Ce = 1, Cg = 0), 
that interacts with a coherent field of the general form 
(10). Substituting the atomic coefficients into (14) 



and 



taking the outer product gives the density matrix 
p of the system. To find the atomic density matrix pa 
we take the partial trace over the field 



Pa 



OO 

n=0 



'n\p\n) . 



(23) 



We can use this result and ( 22 ) to calculate the atomic 
inversion 



(0-3) = Trlpa^s] 



E |CnP cos(2AtV^r+T) 



(24) 



n=0 
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which describes the state of the atom at any given 
time. ((T3) = 1 corresponds to an atom that is fully 
excited, whereas (0-3) = — 1 describes an atom in the 
ground state. Intermediate values imply superposi- 
tions of both states. Figure [T^ shows the inversion 
against scaled time Xt for an excited atom that inter- 
acts with a coherent field with average photon number 
n = 10. The figure clearly shows the collapses and re- 
vivals that are characteristic for the JCM. These are 
often considerered as a proof for the quantised nature 
of radiation [TT] . 

Directly related to the atomic inversion is the expected 
photon number in the field (n), where n = ci^a is 
the number operator. We expect the total excitation 
number N = ^"''^^ + (n) to be constant to conserve 
energy and match the nature of the interaction de- 
scribed above. Direct calculation confirms this ex- 
pectation, giving a graph for (n) that is essentially a 
shifted atomic inversion, mirrored along the (173) = 
axis. 

Another interesting quantity is the von Neumann en- 
tropy of subsystem i 



Si = -Tr[p, log2/5,] 



(25) 



which serves an an entanglement measure for pure 
states |12H16j . giving information about the amount 
of quantum correlations between the atom and the 
field. The von Neumann entropy is less straight for- 
ward to calculate since it involves the logarithm of 
the density operator. To compute this we have to 
diagonalise p, which can be done using Schmidt de- 
composition [5]. The result is shown in Figure [T]d. 
One can see that the system starts with no entangle- 
ment {S = 0) as is expected from a pure initial state 
[17j , but rapidly evolves into an almost fully entangled 
state (5 « 1) and then oscillates at an intermediate 
entropy. Hence, by tuning the interaction time of the 
atom with the field to match exactly the entropy maxi- 
mum, a strongly entangled atom-field system could be 
created, which could then be used as a qubit-pair for 
purposes in quantum information processing (see sec- 
tion V). However, the entropy peaks are very narrow 
and it would prove experimentally difficult to achieve 
the correct interaction time. Note that although we 
have considered the atomic subsystem here, the result 
is generally true for the total system, since it starts 
out in a pure state. Hence the subsystems have equal 
entropy throughout their subsequent evolution [3J |S] . 
For the case of large detuning we found the initial 
state 



\m) 



1 



with Cg = ;^ 



s)) (g) 



a 



(26) 



(for some arbitrary phase 

and the atom in the coherent state to have par- 
ticularly interesting properties. Substituting the co- 
efficients into ( 18 1 gives the state, and hence the full 



find the atomic density operator 
1 r 



Pa 



\9) (.91 + |e) {e\ 

00 

\g) (e|^|C„pe-(^*(^"-i)+« 



(27) 



\e) {g\Y,\C„\'e 



+i(xt(2n-l)+<^ 



It is easy to see from this that the atomic inversion is 
constant, and in fact zero, as is to be expected for the 
off-resonant interaction, since we showed above that 
no excitation exchange takes place between atom and 
field. 

A much more interesting quantity is the entropy. 
Equation ( 27 ) can again be diagonalised by means of 
Schmidt decomposition. The explicit form for the en- 
tropy is then given by 



S{t) = -\g+\' log, |5+p-|5-Plog2l5- 



where 



9± 

Xi 



00 

^|C„|'cos[xt(l-2n)-0)], 

00 

^|C„|'sin[xi(l-2n)-0)], 

n=0 

0. 



(28) 

(29) 

(30) 

(31) 
(32) 



Figure [2] shows the entropy plotted against scaled 
time xt for a coherent state with n = 10. Interest- 
ingly the initially unentangled system rapidly reaches 
a near maximally entangled state in which it remains 
except for interaction times close to xt = nir. This 
result suggests that this particular interaction could 
be used as a reasonably robust method for generating 
entanglement between the atom and the field, since 
it is remarkably insensitive to slight variations in the 
interaction time in the proximity of ~ {n + ^)7r. 
Numerical analysis indicates that the result is inde- 
pendent of the phase </>. 



37! 



FIG. 2: Von Neumann entropy S of an atom in the state 
"75 (Iff) |e)) interacting with a coherent field with n = 
10. 



density operator, at time t. Tracing over the field, we 
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Cavity A Cavity B Detector 



Browne- Plenio Scheme 



Atom 



FIG. 3: Basic experimental setup used for the two 
schemes. An initially prepared atom passes through two 
consecutive optical cavities and is then measured via a 
detector. 



IV. ENTANGLING TWO CAVITY MODES 

In the previous section we described two interac- 
tions that lead to entanglement between an atom and 
an electromagnetic field mode. In this section we shall 
look at a different method that allows for the very pre- 
cise generation of maximally entangled states between 
two optical cavities. 

A. Basic Idea 

The basic concept was initially developed by 
Browne and Plenio [1]. The underlying experimen- 
tal setup is shown in Figure [3j Initially an atom is 
prepared in a certain state. It is then sent through 
two consecutive cavities A and B and allowed to in- 
teract with them for a certain time that is determined 
by the time it takes the atom to traverse the cavities. 
The cavities themselves do not directly interact with 
each other. Once the atom has passed through both 
cavities it hits a detector which makes measurements 
in the {|e) , \g)} basis. This measurement collapses the 
state of the two cavities into a new state, which under 
the correct initial conditions and interaction times, is 
a maximally entangled Bell state. 



B. Generating Bell States 

The four Bell states, given by 
1 



I AB 



I AB 



1 



(|0)^|1)5±|1)^|0)b) (33) 
(|0)^|0)5±|l>^|l>s) (34) 



are the maximally entangled states of a bipartite 
system [H [13 [HI IH [Hj. They are of fundamen- 
tal importance since they form the basis of many 
protocols in quantum information processing, such 
as quantum teleportation which we will discuss in 
section V. Hence it is crucial that we are able to 
efficiently produce these states in the laboratory in 
order to realise these protocols. 



The scheme developed by Browne and Plenio makes 
it possible to generate one specific Bell state, namely 
1^'"'"). In this scheme the atom starts in the excited 
state |e) and both cavities, with frequencies equal to 
the energy sphtting of the atom, i.e. A = 0, are ini- 
tially in the vacuum state |0), giving a combined initial 
state of 



|e)®|0)^®|0) 



B ■ 



(35) 



Since the atom traverses one cavity at a time, first A 
and then B, we can make use of the general solution to 
the time-dependent Schrodinger equation we derived 
earlier in equation (14). First, we can completely ig- 



nore cavity B and only think about the state of the 
system right after the atom leaves cavity A, since we 
know B does not take part in this interaction and 
hence just remains in the state |0)^. Therefore, sub- 
stituting Ce = 1, Cg = 0, Co = 1 and C„ = V n > 



into ( 14 1 we find that the total state of our system is 



\ipinter) = cos(At) |e) |0)^ - i sin(AT) \g) |1)^ 10)^ 

(36) 

after it has interacted with cavity A for a time t = t. 
Now the atom will enter cavity B. Since it now only 
interacts with B and not A, we can redefine our initial 
conditions for this second interaction as 



CL = 



cos(Ar) |0)^ 
-isin(AT) |1)^ 

f 1 if ?i = 
10 if n > 



(37) 
(38) 

(39) 



and use the same method of substitution into ( 14 ) to 
arrive at the solution 



inal/ 



cos2(At)|0)^|0)j 



isin(Ar) cos(At)|0)^|1)s + |1)^|0)b |.g) , 



(40) 



assuming the interaction time is again t = t. 
As described above, we now carry out a measurement 
of the atom in the {|e) , \g)} basis. There are two pos- 
sible outcomes. If we measure the atom in the excited 
state, the cavities collapse back into their initial state 
\tp)cav ~ \^)a\^)b^ which is not a useful state in 
itself, but it allows the scheme to be restarted imme- 
diately. The interesting case is if we measure the atom 
in its ground state. The the two cavity system then 
collapses into the (unnormalised) state 



cos(Ar)|0)^|l)5 + |l)^|0), 



(41) 



By comparing this result to the Bell states in equa- 
tion ( [33| it is obvious that this state approaches the 
state Jl'^) for interaction times Ar — )■ 0. To further 
analyse how close the two-cavity-state is to we 
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FIG. 4: Fidelity F (blue) and success probability P 
(red) plotted against scaled interaction time At for a) the 
Browne-Plenio scheme and b) the scheme we propose. 



we will measure the atom to be in the ground state. 
This probability is 



PbP = 1-\ {e\i^f^nal) |' 

= 1 -cos^(At). 



(43) 



Figure shows the fidelity and success rate against 
scaled interaction time At. The virtue of this scheme 
is its robustness. The fidelity curve is very flat for 
small At, and even at At as high as where Pbp = 
0.75, we still have Fbp ~ 0.97 which could be high 
enough for many applications. However, it can be 
seen that we always have to compromise when using 
this scheme. The fidelity does go to 1 as At — > 0, but 
in this limit the success probability goes to zero. This 
means that we can get arbitrarily close to a Bell state 
by sacrificing success probability, but a perfect Bell 
state is theoretically impossible to achieve. 



2. New Scheme 



can calculate the fidelity F — \ (^caul*^) where we 
first have to normalise our state \iJcav)- 



Fbp — 



1 



cos(At) 



2 cos2(At) + 1 



(42) 



The subscript BP labels this as the fidelity for the 
Browne-Plenio scheme. We see that as expected 
Fbp I as Xt ^ Q. 

However, we also have to consider how likely our ex- 
periment is to succeed, that is, with what probability 



We propose a new scheme, which allows us 
to generate Bell states with arbitrary fidelity and 
non-vanishing success probability. Initially, we 
showed that using the same approach of solving the 
Schrodingcr equation in two separate steps, the cavity 
state after interaction of an arbitrary atom with arbi- 
trary cavity fields of the form of equations ([9| and ( 10 ) 
(the latter of which could be prepared as described by 
Law and Eberly [10]), and after the atom has been 
measured in the ground state, is given by the (unnor- 
malised) state 



oo oo 
n— m— 



71—0 m— 



cosin^^,) sin(f]^) -I- C,C,f+iC^_i M^ti) 
sm{^^) cos(f7^) - CgC^f cosiQ^) cos(f]^)" 



(44) 



B : 



where we have defined 



one described above, with both cavities in the vacuum 



\TA\/n 



(45) 



and an equivalent expression for f2^j. ta and tb 
are the interaction times with cavity A and cavity 
B respectively. It is straightforward to show that on 
substitution of the correct coefficients and by setting 
TA = Tb = T the state (44| reduces to the result ob- 



tained in (41), as expected. 



The new scheme we propose starts with an equally 
simple state as the scheme described above. We sug- 
gest to start with an atom in the ground state (Cg = 1, 
Ce = 0), the first cavity in the number state |1)^ 
and the second cavity again in the vacuum state |0)g. 
This cavity state can be prepared by starting from the 



state as follows. We have derived in equation (36) the 



state of the system after cavity A has interacted with 
an excited atom. If we now measure the atom and 
find it in its ground state, we have prepared the cav- 
ities in the desired state |1)^ lO)^. If we measure |e) 
instead, the cavity reverts back to the original state 
|0)^ and we can try again using another atom. Once 
we have successfully prepared the cavities using these 
excited "auxiliary atoms" we can begin with the ac- 
tual scheme, using a ground state atom. 
The motivation behind the initial state \g) \1) a |0)^ 
is the following. According to the interaction types 
described in section HI, in cavity A the ground state 
atom can cither pass through without any energy and 
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hence excitation exchange with the field, or it can ab- 
sorb the cavity photon, which would bring it into its 
excited state. This leaves the atom in a superposi- 
tion of excited and ground state and the cavity in a 
superposition of |0)^ and Now, as the atom en- 

ters cavity B, its excited component can interact with 
the field by a similar, but in some sense reverse pro- 
cess compared to cavity A, namely deexciting via the 
emission of a photon. This then leaves cavity B also 
in a superposition of |0)^ and 11)^, and in fact the 
combined cavity state can again be turned into a Bell 
state by precise tuning of the interaction time. 
Upon substitution of the coefficients into ( 44 ) we find 



that the cavity will be in the (unnormalisedj state 

\^,av) = sin2(Ar) |0)^ ll)^ - cos'{Xt) 10)^ (46) 

after we measure the atom in its ground state. We 
have again assumed equal interaction times r in both 
cavities. We can compare this to the Bell states in ( 33 1 
to find that \ipcav) = |*~) when sin^(AT) = cos^(AtJ. 
This has the solution Ar — f (2n — 1). The shortest 
interaction time that gives the desired state is thus 
Ar=f. 

Normalising li^cav) we can again compute the fidelity 
of our state in comparison with the Bell state, for an 
arbitrary interaction time. 



F 



1 sin2(2AT) 

2 ^ cos(4At) + 3 



(47) 



Comparing this to ( 42 ) we immediately see one obvi- 
ous advantage of our scheme, namely that the fidelity 
never falls below 0.5. The success probability is given 

by 



P 



1 - ^sin2(2AT). 



(48) 



Both F and P are plotted in Figure ^p. We see 
that although the fidelity maximum and success prob- 
ability minimum coincide, as was the case in the 
Browne-Plenio scheme, our scheme allows the creation 
of states that are arbitrarily close to the desired state. 
The only limiting factor here is the experimental im- 
precision in tuning the interaction time perfectly to its 
ideal value. But assuming that this is a problem that 
can be overcome experimentally, we have a scheme 
that can generate Bell states with a success probabil- 
ity of P = 0.75, as opposed to the vanishing prob- 
ability of the original scheme. The disadvantage of 
our scheme is however the stronger curvature of the 
fidelity curve near its maximum, which implies that it 
is not as robust against variations in the interaction 
time. This is however only an experimental limitation. 



example by cavities of slightly different dimensions or 
by the atom being slightly deflected while travers- 
ing the cavities, hence not taking exactly the same 
straight path in both cavities. Browne and Plenio pro- 
ceed by assuming that the interaction time in cavity 
A is = T and the time in cavity B is = t(1 — e), 
where we have so far not made any explicit assump- 
tion about the size of e other than e < 1. 



Utilising again our general solution (44) we can im- 
mediately incorporate the now different interaction 
times into the previous solution for the Browne-Plenio 
scheme. Thus, finding the new cavity state after the 
atom is measured in its ground state, we now arrive 
at a new fidelity 



P'bp 



1 



cos(At) sin(Ar) sin(AT(l — e)) 



(49) 



2 cos2(At) sin2(Ar(l - e)) + sm^iXr) 
A similar treatment of our scheme gives the fidelity 



F' = 



1 sin(2AT)sin(2AT(l - e)) 

2 ^ 2cos(2AT)cos(2Ar(l - e)) -i 



(50) 



The two fidelities are plotted as functions of At and 
e in Figures [5^ and [5]3 respectively. We can see 
that both schemes exhibit extended areas for which 
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C. Interaction Times 

Browne and Plenio also considered the effect of dif- 
ferent interaction times in the two cavities due to ex- 
perimental imprecision [T]. This could be caused for 



FIG. 5; Contour plots of the fidelity as a function of inter- 
action time and deviation e for a) the Browne- Plenio 
scheme and b) the scheme we propose. The red contour 
corresponds to a fidelity of 0.99. The remaining fidefity 
contours are at 0.1, 0.2, 0.9. 
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F > 0.99, which implies a certain robustness against 
experimental inaccuracies. In line with Figure [4| we 
see that the Browne-Plenio scheme has a flatter fi- 
delity function. The contour plots also demonstrate 
the interesting possibility to use different interaction 
times for cavity A and B to achieve the desired states. 
In fact our scheme allows us to find a suitable e for 
any interaction time ^ < tt. 

D. Thermal States 

We can further extend both schemes by asking what 
would happen if we allowed the system to thermally 
interact with its environment. This can be done by 
allowing either the atom or the cavities to be in the 
thermal state, which is described by the density oper- 
ator 



PTh 



(51) 



where /? = {kgT) ^ with kg being the Boltzmann 
constant, T the temperature, and H the atomic or 
field Hamiltonian discussed in section II [TU]. For the 
atom it can be shown jSl [5] that this is equal to 



PTh 



with 



= |e) (e| + \g) {g\ , 
Pg = (e-^""+l)-i. 



(52) 

(53) 
(54) 



A similar result can be obtained for the thermal states 
of the field, where 



PTh 



with 



and 



Pn = 



(1 + n)"+i 



(e*^" - 1)- 



(55) 



(56) 



(57) 



Both these thermal states are distinctly different from 
the states we have considered so far. They are intrin- 
sically mixed states, representing statistical mixtures 
of all the accessible states. Hence they can not be sep- 
arated into pure states as was the case with the initial 
states considered so far. The only way to describe the 
thermal states is via the density operators introduced 
here. 

This means we can now no longer rely on the solu- 
tion of the Schrodingier equation (14 1, but have to 
make use of the evolution operator ( |16[ ). The density 
operator then evolves according to [10] 



p{t)^U{t)p{Q)ij\t). 



(58) 



1. Browne-Plenio Scheme 

First we look at the Browne-Plenio scheme where 
we now allow the two cavity modes to be in thermal 
states. The atom still starts out in the pure state 
|e) which can be described by the density operator 
Pa = |e) (e|. We can again treat the two interactions 
individually. The density operator for the joint system 
of the atom and cavity A is 



/5(0) = /3„ (g) pI^^^ 




(59) 



where the second line refers to the atomic basis. We 



can now use equation ( 58 ) to find the density operator 



of the system after the atom has interacted with the 
thermal field in cavity A for t = r. The resulting 
4 X 4-matrix can be defined as 



P{r) = 



Peg 

Pgg 



(60) 



The total density matrix before the interaction with 
cavity B is then p{t) ® p^^^ , i.e. 



p'(0) = [P<'<'® PThs P<^9 ® PThE 
\Pge®PThs Pgg® PTh, 



(61) 



in the atomic basis. Using the evolution operator to 
once more evolve this state by a time r gives the final 
state of the full system, p final- When we measure 
the atom in the ground state, the remaining cavity 
density matrix is p^av = {9\Pf^nai\g)- The result can 
be written as 



= {Cp^^C)a ® {Sp^J')b 
+ {Sp^hC)A ® (C'p^J')b 
+ {Cp^J% {Sp^,,C')b 
+{Sp^hS')A ® {C'p^hC')i 



(62) 



where C, S, C and S" are the matrix elements of the 
evolution operator U — 



c s' 
s c" 



(c.f. equation ( 16 1). 



The exact evaluation of pcav following the procedure 
above is straightforward albeit very laborious and is 
left out here. 

After truncating the series describing the thermal 
states for n > 1 (since then the P„ are of negli- 
gible size under normal conditions, see below), pcav 
can be expressed as a 9 x 9-matrix in the field basis 
{|0)^|0)5,...,|2)^|2)5}. 

At this point we are no longer able to simply compare 
this state with the Bell state to judge how much en- 
tanglement we have. Furthermore, the von Neumann 
entropy is not a well defined quantity anymore since 
we are dealing with mixed states [5]. We have to rely 
on more sophisticated entanglement measures. One of 
the simplest such entanglement measures is the loga- 
rithmic negativity Em{p) [131 [HI [12] ■ It is based on 
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a) ' 




b) ' 



higher than this does i?jv substantially fall off. Hence 
if we assume that our cavity operates in the microwave 
range of the spectrum with luq « 10 GHz, we can op- 
erate our cavities at temperatures of up to T « lOOK 
and still achieve very high levels of entanglement. We 
can also use these values to justify truncating the ther- 
mal states for n > 1. For the highest temperatures 



shown in Figure 
ing to (57) n w 
Po « 0.53 > Pi ? 



^ = 10 



K 
THz ■ 



we have accord- 



87. Substitution into (56) gives 
0.25 > P2 



0.11. For the region 



we are particularly interested in, 



< 1 



K 



these 



Scaled Tempeiatiire ^ {KfTHz) 



T_ 

ujQ ^ ^THz 

results are even more dramatic, with n 4.8 x 10^^ 
and 

Po « 0.9995 > Pi w 4.8 x 10"'' > P2 w 2.3 x 10"^ 

(66) 

which justifies the approximation made by omitting 
terms containing P2 or higher. 



2. Our Scheme 



FIG. 6: Logarithmic negativity En as a function of the 
scaled temperature ^ for a) the Browne- Plenio scheme 
with an interaction time Ar = O.OItt and b) the scheme we 
propose with Ar = 0.2887r. 



the negativity 



lA.I-A, 



(63) 



where the Xi are the eigenvalues of the partial trans- 
pose of p = ^Pijki K) (j| ^ \k) {l\ over subsystem A 
defined as 



i,j.k,l 



(E>\k) {l\ 



(64) 



The logarithmic negativity can be calculated from the 
negativity by 



(65) 



It can be shown that the result does not depend on the 
subsystem over which the transpose is taken [12l HI] . 
Its value falls in the range < Ejsi < 1 where is a 
pure state and 1 a maximally entangled state. Ei\[ is 
also an upper bound to the distillable entanglement 
Ejj, a result to which we shall come back in our later 
discussion of teleportation. 

After taking the partial transpose over cavity A, we 
numerically calculated the logarithmic negativity of 
the state Pcav Figure |6^ shows E]\[ after an inter- 
action time At = O.OItt, i.e. very close to the ideal 
time At — >■ 0, plotted against the scaled temperature 
^ . For low temperatures ^ ^ 1 yth logarithmic 
negativity is very flat and almost equal to 1 imply- 
ing near-perfect entanglement. Only for temperatures 



For the scheme we proposed, we analysed the case 
where the atom is thermally interacting with its envi- 
ronment and the cavities are in the state |1)^ \0)g as 
before, i.e. p(0) = p^;^ |1) (1|^ ® |0) (0]^ • Since 
only contains two terms, we do not have to make any 
approximations as in the previous case. Following the 
same procedure outlined above for the Browne-Plenio 
scheme, we arrive at the density matrix for the cavities 
after the atom has been measured in \g) 

o\ 






0/ 

(67) 

in the basis {|00) , |01) , |10) , |11) , |20) , |21)}, where 
we have introduced the shorthand notation 




















PA 













-PgClsj 


PA 



















PeC2S2Sl 











PeC2S2Si 


P s2 


Vo 















c„ = cos(ATvn), 
s„ = sin(AT\/n). 



(68) 
(69) 



This can now in fact be directly compared with the 
Bell state density matrix p^- = 1^ ) I in the basis 
{|00),|01),|10),|11)}. 



/O o\ 

+i -i 

-| +1 

\0 0/ 



(70) 



We find that there are two conditions that have to be 
satisfied for our scheme to produce a Bell state. The 
first condition is Pe — >■ 0, i.e. T — > 0, which was to 
be expected since the ideal scheme assumes that the 
atom is in the ground state. The second condition is 

sin'^(AT) = cos2(At) = cos(At) sin2(AT) (71) 
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which has the real solutions 



At = 2 ( nn ± tan 



7^-2), 



(72) 



with the lowest positive solution Ar w 0.2887r, which 
is slightly longer than the ideal value of the thermally 
isolated scheme, Xt — j. 

Further, we calculated the logarithmic negativity for 
this scheme in the way described above. The results 
are plotted m^p for an interaction time At = G.2887r. 
The results look remarkably similar to the results 
found for the Browne-Plenio scheme. The low temper- 
ature behaviour is equivalent for both, having a near 
perfect entanglement for ^ < ^Yff^- O^^'' scheme 
appears to be considerably less sensitive to increas- 
ing temperature and hence more reliable. However, 
a direct comparison is difficult since in one case the 
field modes and in the other the atom are allowed to 
interact thermally. 



V. APPLICATIONS TO TELEPORTATION 

Finally, we want to analyse the possibility of apply- 
ing our scheme to achieve quantum teleportation |23j . 
Let us assume that once we have created the entan- 
gled state pcav we can transmit the two cavity modes, 
say cavity mode A to Alice and B to Bob, without the 
state decohering. We now want to make use of the en- 
tanglement that Alice and Bob share to transmit an 
unknown qubit |V-')(7 



|V)c = a|0)c 



611) 



(73) 



that Alice possesses to Bob. We want to transmit the 
state using only local operations and classical com- 
munication (LOCC). This is the concept of quantum 
teleportation [51P^[Tgj. 

Let us first assume that Alice and Bob share a perfect 
Bell state 1^*")^^ ^^"^ leave the generalisation to our 
imperfect state for later. The total three-qubit system 
is in the state 



/AB 



(74) 



By making full use of the Bell states ( 33 ) and ( |34[ ) we 
can rewrite this state as 



Ix)-^[|*+)ac®("«|0>b + H1)b) 



'AC 



6|1) 



|$+)^P® (-fa 11)^-6 10)5) 
I$")ac®(+«|1)b+M0)b) 



(75) 



collapse to one of four possible states, depending on 
the outcome of Alice's measurement. The results are 



I*-; 
1$-; 
1$+; 



AC 
AC 
AC 
AC 



-«|0)b 

-all 



-a 



b\l] 
b\0] 



-a|l)„-6|0) 
611! 



B ' 
B ' 
B ' 
B ' 



(76) 
(77) 
(78) 
(79) 



where the notation |a;) — > \y) stands for Alice measur- 
ing |a;) and Bob's state collapsing to \y). Now Alice 
can tell Bob via a classical communication channel 
what she measured, which will tell Bob the state of 
his qubit. If Alice measured qubit immedi- 

ately takes the form |x), thus completing the telepor- 
tation process whithout Bob having to do anything. 
However, if Alice measures the states 1$^)^,^, 1*1'^)^ 



/AC 



/AC 



, Bob can use the local operations a^, i&y 
or cTz respectively (the ct^ being the usual Pauli matri- 
ces) to transform his qubit into the desired state \x), 
thereby also completing the teleportation process [H] . 
In the process, the entanglement that Alice and Bob 
shared is "used up". There are no correlations be- 
tween Bob's qubit B and Alice's two qubits A or C 
anymore. However, A and C are now in a Bell state, 
but they are only locally available to Alice, so they 
cannot be used for further teleportation between Alice 
and Bob. It thus makes sense to speak of entangle- 
ment as a resource. 

Now that we have worked out the basic teleportation 
protocol we can analyse what would happen if Alice 
and Bob do not share a perfect Bell state but the 
state Pcav given by equation (67 1. The total state of 



the three-qubit system is described by the density op- 
erator 



Ptot Pcav 



(80) 



Alice now makes the same measurement in the Bell 
basis on her two qubits A an C as before. She thereby 
projects Bob's qubit into one of the following states 



Pb 
p'b 
Pb 

pi 



= {■^~Ac\Ptot\^~Ac) ^ 

= {^-Ac\Ptot\^~Ac) ^ 

- {^\c\Ptot\^Ac) ^ 

= {-^Ac\Ptot\^Ac)- 



(81) 

(82) 
(83) 
(84) 



Bob then has to apply the same local transformations 
given by the Pauli matrices as in the ideal case de- 
scribed above. The fidelities of the final teleported 
states compared to the original state, after Bob has 
made the necessary transformation, are given by 



So far there was no action by Alice or Bob involved, 
only a change of basis in which we express our state. 
Now we can immediately see what happens when Alice 
makes a local measurement on her two-qubit state AC 
in the Bell basis. The state that Bob possesses will 



Fa 

Fb 
Fc 
F, 



Io'zPb'^z 



(85) 
(86) 
(87) 
(88) 
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FIG. 7: Fidelities Fa 
function of scaled temperature 



Fd of the teleported state as a 
3_ (in ^) and \a\\ The 
interaction time is At = 0.2887r, close to the ideal value. 



We calculated the explicit forms of these fidelities and 
found 



'4 



-2|a|2|6|2cis2 + |5|4,4)^p^|^|2|^| 



12^2 



2„4 



(89) 



and 



Pg{\a\^si + 2\am\,sl + \b\^ci) + PMW4sl 



|a|2s4 + |fo|2c2 



(90) 



We can eliminate one variable by insisting that the 
state \ip) is normalised, implying |a|2 + |6|2 — l. The 
fidelities are plotted as functions of temperature and 
|a|2 in Figure [t] for the near-ideal interaction time 
At = G.2887r. The plot explicitly shows Fa ~ Fd, but 
the other plot for Fi, — Fc is essentially indistinguish- 
able and was therefore omitted. In fact, a numerical 
analysis shows that Fa — Ff, — Fc — F^ for At going 



to one of the values in equation ( 72 ) 



Surprisingly, we see that now the fidelity of the tele- 
ported qubit has a direct dependence on the precise 
nature of the original qubit, through the |a|2 depen- 



dence in equations (89) and (891, unlike in the ideal 



condition protocol, which is completely independent 
of the explicit form of \ip)- Yet, the figure shows that 
this effect is barely noticeable and converges to the 
correct value F = 1 as T — )■ for all |a|2. We also 
note that the fidelity drops with increasing temper- 
ature, in line with the results found for the thermal 
states in section IV. D. In fact, the temperature depen- 
dence of the fidelity shows a remarkable resemblance 
to the temperature dependence of the logarithmic neg- 
ativity shown in Figure |6] As was to be expected, if 
we want to achieve perfect teleportation using only 



one entangled qubit pair produced by our scheme, we 
have to go to infinitesimally low temperatures. 
However, there exist more realistic ways of realising a 
perfect teleportation protocol. It can be shown [T31,!T^ 
that given a certain number N of non-maximally en- 
tangled qubit pairs such as our state (67), it is possible 



to "purify" the entanglement, producing M < N max- 
imally entangled states, using local operations only. 
The crucial quantity limiting the efficiency of this pu- 
rification process is the distillable entanglement Ed 
briefly mentioned in the discussion of the logarith- 
mic negativity. The condition limiting the efficiency 
is M < NEo- We can turn this into an inequality in- 
volving the logarithmic negativity calculated earlier, 
since as mentioned above the logarithmic negativity 
is an upper bound on the distillable entanglement. 
En > Ejj. This implies 



N > 



M 



ENipy 



(91) 



Hence if we want to use our scheme to provide M 
maximally entangled qubit-pairs that can then each 
be used to teleport one qubit, we need at least 
Nm.iri = p fi — T copies of the state pcav shared 



between Alice and Bob. Note that this only allows us 
to determine a lower limit on the necessary number 
of copies. There is no guarantee that this number will 
be sufficient to actually achieve the desired number of 
maximally entangled states through the purification 
process. Nevertheless, these considerations at least 
qualitatively show that it is possible to use our 
scheme to generate entangled states that can be 
sent to two parties Alice and Bob, which can carry 
out entanglement purification procedures, and then 
use the resulting states for teleportation with unit 
fidelity. Considering again the results in Figure [6] we 
can see that even at relatively high temperatures it 
is possible to provide states for the teleportation pro- 
tocol, as long as we are willing to sacrifice efficiency. 
Quantitatively we can only place a lower bound on 
the required number of copies needed. Under real 
conditions this is likely to be increased even further 
by the fact that the states could undergo further 
decoherence as they are distributed to Alice and Bob. 



VI. CONCLUSIONS 

We have derived the Rabi model Hamiltonian and 
justified its simplification under the rotating wave ap- 
proximation into the Jaynes-Cummings model Hamil- 
tonian. These results lead us to show that the JCM 
predicts two distinct interaction processes between the 
atom and the field. Near resonance the two subsys- 
tems can exchange energy, whereas for large detuning 
only the phases of the two subsystems are affected. In 
section III.B we showed how both of these could be 
used to entangle an atom-field system. 
We then extended the scheme developed by Browne 



12 



and Plenio for entangling two cavities by proposing 
a slightly modified scheme, promising a considerably 
higher success rate and a fidelity only limited by the 
(im)precision of the experimental setup. 
We further added realism to both schemes by allow- 
ing the systems to thermally interact with the envi- 
ronment and showed that in this case the fidelities are 
monotonously decreasing functions of temperature. 
Finally we considered the applicability of our scheme 
to generating states that can be used for quantum 
teleportation. We showed that the fidelity of the tele- 
ported state also monotonously decreases with tem- 
perature. Nevertheless, perfect teleportation is pos- 
sible if we add a purification scheme to the protocol, 
which would in turn require a larger number of shared 
states to be generated and distributed. We were able 
to place a lower bound on this number related to the 
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